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We study a one-dimensional chain of harmonically coupled units in an asymmetric anharmonic 
soft potential. Due to nonlinear localisation of energy, this system exhibits extreme events in the 
sense that individual elements of the chain show very large excitations. A detailed statistical analysis 
of extremes in this system reveals some unexpected properties, e.g., a pronounced pattern in the 
inter event interval statistics. We relate these statistical properties to underlying system dynamics, 
and notice that often when extreme events occur the system dynamics adopts (at least locally) 
an oscillatory behaviour, resulting in, for example, a quick succession of such events. The model 
therefore might serve as a paradigmatic model for the study of the interplay of nonlinearity, energy 
transport, and extreme events. 


I. INTRODUCTION 

Extreme events as naturally occurring large deviations of 
a system’s state from its normal behaviour have gained 
considerable interest during the past years [1-4]. Many 
phenomena where one likes to speak about extreme 
events are related to natural disasters, such as earth¬ 
quakes, weather extremes, rogue waves. Due to the pos¬ 
sible consequences of such events, much work has gone 
into the prediction of where and when they will occur. 
This allows for better planning in the affected regions; 
for example, better evacuation strategies in coastal areas 
where tsunamis are a threat. 

However, also in physical models extremes have been 
discussed, e.g. in microwave cavities [5], optical fibres 
[6], superfluids [7], and semiconductor lasers [8], to name 
but a few. When studying extremes, four issues are in 
the main focus: the statistics of extreme events, in par¬ 
ticular the question of how frequently very large events 
are to be expected; temporal correlations between succes¬ 
sive events; the underlying mechanism for the creation of 
such large deviations; and the issue of how to forecast 
individual events. 

An extreme event in a spatially extended system often 
presupposes that a significant amount of energy localises 
in a small area of a system (with respect to the systems 
overall size). This problem of energy localisation has re¬ 
ceived much attention in the area of discrete lattices. The 
genesis of this focus comes from work carried out on what 
is now referred to as the Fermi-Pasta-Ulam problem [9]. 
One of the remarkable discoveries from their numerical 
work on a one-dimensional lattice system consisting of 
nonlinearly coupled units, was that discrete systems with 
added nonlinearity can localise energy. This led to excit¬ 
ing work on discrete breathers (also known as intrinsic 
localised modes) [10, 11[, which are large amplitude, ex¬ 
ponentially localised, exact periodic solutions to nonlin¬ 
ear discrete lattice systems. Such solutions are long-lived. 
In contrast, it is also possible for these lattice systems to 
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exhibit transient dynamics which have the characteristic 
of exponentially localised energy; for example, chaotic 
breathers [12, 13[. These solutions typically have finite 
lifetime with the energy eventually dispersing to other 
parts of the system. 

The so-called extreme events have been observed in 
a number systems; for example, in a one- and in a 
two-dimensional discrete nonlinear Schrodinger system 
[14, 15]. The exact mechanism(s) generating extreme 
events are difficult to pin-point, and thus these stud¬ 
ies often take the form of a statistical analysis. How¬ 
ever, one mechanism that is known to account for energy 
localisation, and indeed a likely candidate as a mecha¬ 
nism that produces extreme events, is known as stochas¬ 
tic localisation. This mechanism is particularly present 
in systems where the oscillators under investigation are 
soft, i.e. those whose frequency of oscillation decreases, 
and the density of states increases, as their amplitude in¬ 
creases. Thus, it is favourable (entropically speaking) for 
such oscillators to seek out regions of phase-space with 
the larger density of states [16, 17[. 

In this paper we study a model system where we high¬ 
light a mechanism which might be typical of extremes: 
the localisation of energy at some few degrees of freedom. 
We present a detailed statistical analysis of extremes in 
this system with some unexpected properties, e.g., a pro¬ 
nounced pattern in the inter event interval statistics. The 
model being Hamiltonian, the total energy is conserved. 
This might be relevant even for real world phenomena, 
i.e. in dissipative systems, such as extreme precipitation 
events: also there, mass conservation and energy con¬ 
servation play a role, since there cannot be more rain 
than water vapour transported by air, and latent energy 
has to be released when water vapour condensates into 
rain drops. We also examine the differences and similar¬ 
ities in the dynamics and statistical distributions when 
using two different thresholds that determine the onset 
of extreme events - one threshold being an amplitude 
Q = (Ithresi ^nd the other being a local energy (potential 
plus kinetic) threshold E = Ethres- Further, we explore 
how the system reacts to changes in the coupling strength 

K. 
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II. SYSTEM 


The system under investigation is a one dimensional con¬ 
servative system consisting of N coupled units which has 
a Hamiltonian given by 


N 


H{p,q} = 


n—1 


1 Kj 

^pI + U{qn) + ^{qn+i - qnf 


( 1 ) 


where Pn and represent canonically conjugated mo¬ 
mentum and positions variables, and with local onsite 
potential 


U{qn) = D{l-e-^‘>’'f ( 2 ) 

which is known as the Morse potential [18]. The parame¬ 
ters D and a control the potential well height and width 
respectively. In the limit g —>■ oo, U{q) D. Taking the 
following dimensionless variables q'^ = ag„, t' = ay/Dt, 
and k' = K,/a?D, the corresponding rescaled equations of 
motion are given by 


g„ = -2(l-e '^")e «"-|-K(g„+i-|-g„_i - 2g„) (3) 

where we have dispensed of the prime notation. This 
system now has a single parameter k which regulates 
the strength of the coupling between neighbouring units. 
Throughout we will be using periodic boundary condi¬ 
tions and, unless otherwise stated, a chain will consist 
of 100 units. Moreover, the system’s energy density will 
be equal to unity, i.e. '^En = N. In general terms, 
the anharmonic onsite potential can be seen to repre¬ 
sent intramolecular interactions, while the linear cou¬ 
pling between the units represents intermolecular inter¬ 
actions. The particular onsite potential considered here, 
the Morse potential [18], is regarded as a suitable model 
for the potential energy of a diatomic molecule because it 
can account for the anharmonicity of real chemical bonds, 
and it implicitly includes the effects of dissociation (bond 
breaking). 

To begin, there are two limits that should be explored. 
First, the uncoupled limit k —^ 0 will exclude the possi¬ 
bility of energy propagating through the chain, with each 
chain unit maintaining its initial energy. Thus, if some 
chain unit possesses energy > 1 then this unit can 
travel to infinity in the positive q direction. The other 
limit, K —> oo, results in infinitely strong interactions be¬ 
tween neighbouring units. For systems with finite energy, 
the only possible chain configuration is a uniform rod-like 
state where the dynamics of each unit is the same. Any 
deviation from this rod-like state would put an infinite 
amount of energy into the harmonic springs, something 
that is clearly not possible in the finite energy case. Im¬ 
portantly, the limits k = 0 and k = oo are integrable. 
The interesting collective dynamics, and consequently 


the interesting statistics, come in the nonintegrable case, 
i.e., for finite nonzero values of the coupling strength. 

In this paper we are primarily concerned with extreme 
events. An extreme event is said to have occurred at 
some chain location (unit) when this unit has an ampli¬ 
tude greater than a predefined threshold, that is when 
qn > qthres = 4. Here, we will also look at an energetic 
threshold defined in terms of local energy. An extreme 
event in terms of local energy occurs when the sum of the 
potential plus the kinetic energy of a chain unit crosses 
a threshold, that is i5„ > Ethres = 3. The coupling part 
of the potential energy is shared evenly between neigh¬ 
bouring units. 

We expect that small changes in the above thresholds 
will not have qualitative implications for the statistical 
results that will be discussed later. Of course, other 
thresholds will naturally produce different results, and 
there is no perfect way to define the thresholds. In¬ 
deed, we are fully aware that the statistics discussed 
later depend on how the thresholds are defined. How¬ 
ever, the thresholds chosen here seem reasonable. The 
g = 4 threshold, being well beyond the inflection point of 
the (Morse) potential at g Ri 0.693, corresponds to this 
potential having a gradient close to zero, with energy 
(excluding the coupling and kinetic energies) E fa 1 = 
U{oo). The E = 3 threshold means that a unit will have 
energy that is at least 3 times the energy density of the 
system. Moreover, choosing a fixed threshold allows for 
a direct comparison between the statistics for different 
coupling strengths. Thus, an extreme event occurs when 
a unit either crosses a particular point in space, or when 
it attains a sufficient amount of energy. These thresholds 
are quite different, and importantly, a threshold crossing 
in space or energy does not necessarily imply a crossing 
in the other; that is, a high concentration of energy at a 
particular chain unit, does not imply a large amplitude 
excitation, and vice versa. 


III. DYNAMICS 

With a view to the later discussion on extreme events 
it is worthwhile to look at the dynamics of the system un¬ 
der investigation. For this, we inject kinetic and potential 
energy randomly into the system such that £’„(0) = 1 for 
all n, and subsequently study the evolution of this initial 
excitation. Fig. 1 contains heat maps which allow one to 
observe the evolution of the amplitudes g„ through time. 
Interestingly, it can be seen that chain units can undergo 
large fluctuations in their amplitudes. In the correspond¬ 
ing energy heat maps (not shown) it can be seen that the 
energy can become concentrated (localised) at certain 
chain locations. The localised island-like structures, in 
amplitude and energy, that emerge are able to move along 
the chain; the speed of this mobility seems to depend on 
the strength of the coupling between neighbouring units, 
the number of units involved in the structure, and on the 
energy contained in the structure. This mobility allows 
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Figure 1. Color online. Heat maps showing the evolution of 
the amplitudes qn{i) for 100 coupled units, as a function of 
time, for two different values of the coupling strength: n = 
0.1 (top), K = 0.4 (bottom). All amplitudes q-nit) > 4 are 
coloured blue. Here t denotes time, and n the location of 
the units on the chain. The oscillating island structures are 
clearly visible. 


for the interaction of two or more localised structures re¬ 
sulting in various consequences. A number of possible 
outcomes from these interactions include, the creation of 
a larger structure through the merging of smaller ones, a 
dispersion of the energy contained in the structures, two 
structures reflect off one another, etc [19]. This presents 
an algorithmic challenge when trying to monitor the time 
evolution of these localised structures, as it is unclear 
what happens to the individual structures when they col¬ 
lide with others. Looking more closely at the individual 
localised structures, it can be seen that they are oscillat¬ 
ing in time. This can be observed in Fig. 1 by noticing 
that the localised structures with high amplitudes are in¬ 
termittently interrupted by periods of low amplitude, i.e. 
in the heat maps the colours of the localised structures 
alternate between green and blue. 

The system dynamics also react in interesting ways 
when the coupling strength is changed. With increasing 
coupling strength energy can disperse more freely along 
the chain. Tims, one can see localised structures mov¬ 
ing over a significant portion of the chain in short pe¬ 
riods of time. On the other hand, increasing the con- 
pling impedes the formation of large amplitnde localised 
structures. This is because more energy is reqnired in 
the harmonic springs. This, in turn, means that there 
will be fewer qthres threshold crossings with increasing 



Figure 2. Number of qthres (filled circles) and Ethres (crosses) 
threshold crossings as a function of the coupling strength k. 


coupling strength. Fig. 2 shows the number of qthres 
threshold crossings for the chain units as a function of the 
coupling strength k for an ensemble of 1000 initial condi¬ 
tions, and simulation time per initial condition t = 5 x 10^ 
(si 1.1 X 10^ X To with To = being the period du¬ 
ration of harmonic oscillations about the potential mini- 
mnm). The cnrve representing the nnmber of threshold 
crossings shows a tendency, as k is increased, to approach 
a constant (see below). However, one needs to consider 
that in the uncoupled case k = 0, for the initial condi¬ 
tions considered here, each unit in the chain will cross the 
threshold once before escaping to infinity. This results in 
100 X 1000 = 10® (number of units times number of initial 
conditions) threshold crossings. Therefore, moving into 
the conpled regime, one must conclude that there is a 
transition which sees a dramatic increase in the number 
of extreme events in q before this behaviour reverses and 
there is a decay in the number of events. One may be 
inclined to think that there is a singularity at k = 0 caus¬ 
ing a large jump in the number of events as one moves 
away from this limit. However, the finite simulation time, 
together with the very weak coupling that allows chain 
units to travel large distances almost unhindered, more 
likely points to a gradual smooth increase in the number 
of events as k is increased from k = 0. 

Dynamical behaviour, analogous to that discussed 
above, can be observed for energy heat maps (not shown), 
where the local kinetic plus potential energy for each 
chain nnit are monitored. However, the oscillating na¬ 
ture of the localised structures is less apparent. This is 
due to the fact that energy in a unit can stay the same 
even when its amplitude drops, as the potential energy 
can be converted into kinetic energy. The consequences 
for the number of Ethres threshold crossings as k is in¬ 
creased can be seen in Fig. 2. Notably, the number of 
crossings increases with the coupling strength. Again, it 
is easy to deduce the number of crossings in the k = 0 
case: without coupling to other units, no unit is able to 
attain the additional required energy to cross the thresh¬ 
old, and so there will be no energetic extreme events. 
Further, while large amplitude excursions by the chain 
units become less likely with increased k, the same can¬ 
not necessarily be said for the local energy. That is, local- 
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isation of energy is likely to be enhanced with increased 
coupling strength. Of course this is true only up to some 
critical k value, as the chain will approach the rod-like 
configuration (mentioned earlier) due to energetic con¬ 
straints for high values of k. In fact, this critical value is 
almost certainly k = oo (discussed below), as for any fi¬ 
nite coupling strength, there is nothing that forbids small 
amplitude deviations between neighbouring units. This 
presents a clear difference between the two chosen thresh¬ 
olds. That this same process (of increasing k) allows for 
more frequent energetic extreme events, in contrast to 
the qthres extreme events, is interesting. It would seem 
that energy still localises efficiently without the need for 
the chain units to attain large amplitudes. 

Finally, it is worth discussing the number of qthres and 
Ethres threshold crossings in the case k = oo. We re¬ 
iterate, that in this case, the only possible chain con¬ 
figuration is that where the chain is rod-like and each 
unit undergoes exactly the same dynamics. This means 
that zero energy is contained in the coupling part of the 
potential. Now we have a chain acting (effectively) like 
a single unit, with energy E = 1. This is exactly the 
K = 0 case, and so the number of threshold crossings, 
qthres s-iid Ethres, are exactly the same in both the k = 0 
and K = oo cases: = n(g^)(;“) = 10^ and 

= n{E^f;;;.^^) = 0. Turning our attention to the 
curves contained in Fig. 2, the g-curve will likely con¬ 
tinue with an exponentially slow decay, with increasing 
At, towards its At = 0 value. On the other hand, the E- 
curve must at some point reverse its behaviour and begin 
a descent towards zero with increasing At. 


IV. ENERGY LOCALISATION 

To give a quantitative characterisation of how well the 
energy localises for various coupling strengths, we intro¬ 
duce the following quantity [12]: 


LN{t) 



N 




Efjt) \ 

m)f/ 


( 4 ) 


where N is the number of units in the chain, and Ei is 
the local potential plus kinetic energy contained at chain 
location i. The angled brackets indicate an ensemble 
average over 100 initial conditions. For full equipartition 
of energy, that is Ei = E/N for all i and total system 
energy E, Lm = 1. Conversely, if all the system energy 
is concentrated at some single chain location n, that is 
En = E and Ei = Q for all i ^ n, then = N. 

The results for various values of the coupling strength 
K are contained in Fig. 3. Starting from a state of full 
equipartition of energy a clear trend emerges, regardless 
of the coupling strength. The early evolution of is 
characterised by an overall increase from its initial value 
of E]sj = I. Subsequently, Ljv settles into its long term 
behaviour plateauing close to a maximal value. It should 
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Figure 3. Color online. The curves at the top of the figure 
show the evolution of the localisation measure Ln, the middle 
curves are the potential energy, and the bottom curves are 
the kinetic energy. For each coupling strength (see key) the 
curves have been obtained using the same ensemble of 100 
trajectories. Note that there is a break in the axis at t = 45. 


be noted that the evolution of Ljv for individual trajec¬ 
tories is much more erratic with numerous peaks which 
go significantly above the average values. Also worth re¬ 
membering is that Ljsi remains constant in the case of 
uncoupled units - that means Ljv = 1 for the initial set 
up considered here. While the measure of localisation 
might appear low, L]\[ N, even those events which 

show large amplitude excitations need to be considered 
with respect to a high energy background. This will be 
discussed in the coming paragraphs. 

The evolution of Lav can be related to the correspond¬ 
ing evolution of the average local kinetic and potential 
energies (see also Fig. 3). In particular, a transient pe¬ 
riod, where the average kinetic and potential energies os¬ 
cillate erratically, sees an increase in Ljv, i.e. an increased 
localisation of energy. Afterwards, the soft nature of the 
underlying Morse potential emerges with the average po¬ 
tential energy becoming higher than the corresponding 
kinetic energy. Moreover, these energies now undergo 
small oscillations about some constant value which is de¬ 
pendent on the coupling strength. At the same time, L^ 
has settled into his long term behaviour. 

To see how the energy localises for extreme events we 
show in Fig. 4 an excitation pattern (with the positions 
centred at n = 0) obtained by averaging over the max¬ 
imum amplitude of all energetic extreme events, and its 
30 nearest neighbours, for a collection of 2 • 10® extreme 
events. To be precise, let im(t') be the position of the 
largest excitation at the time t' of the extreme event, 
with Ei^{t') > Ei{t') Vi ^ im'. we find the mean of 
2 • 10® such excitations with the position shifted to 

zero. Whereas the linear scale shows a strong concentra¬ 
tion of the excitation pattern at 0, which is the position 
of the maximum, a logarithmic representation, after sub¬ 
tracting the background E^i^ from each curve, shows ev¬ 
idence for exponential localisation. Qualitatively similar 
results have been obtained for the averaged g-amplitude 
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Figure 4. Color online. The average energetic excitation pat¬ 
tern in the vicinity of the location of the maximum im, when 
the excitation exceeds the Ethres = 3 threshold. The logarith¬ 
mic representation (inset) shows that it is, to a good approx¬ 
imation, the superposition of some constant background and 
a symmetric exponential part, hence we observe exponential 
localisation in space. 


excitation pattern (not shown). 

Returning to the localisation measure L^r, we find 
that the curves shown in Fig. 4 fall in the range Rat S 
(1.2,1.35). However, when the constant background 
Emin is removed from each curve, L^r then falls in the 
range € (12.8,19.2) - that is 1 <C Rat < N = 31. 
This demonstrates that all energy localisation must be 
considered in terms of a high energy background, and 
thus explains why the curves for Lat in Fig. 3 appear to 
suggest that energy does not localise efficiently in this 
system. 


V. STATISTICS OF THRESHOLD CROSSINGS 


To gain further understanding of threshold crossings, 
going beyond an analysis of the internal dynamics which 
produce such crossings, it is useful to examine their sta¬ 
tistical properties. We consider an event to be extreme 
when a chain unit crosses over the threshold qthres = 4 
(Ethres = 3). This event begins with a threshold crossing, 
i.e. qn(ti) > qthres-, and ends when this same units re¬ 
crosses the threshold, i.e. qn(t 2 ) < qthres, where ti and t 2 
denote the times of the threshold crossings (and similarly 
for En). Fig. 5 contains a number of distributions detail¬ 
ing various statistical aspects of threshold crossings, for 
and En respectively, and in particular how these dis¬ 
tributions change as a function of the coupling strength 

K. 

For the simulations we used 1000 initial conditions, 
each of which was chosen randomly, but with the restric¬ 
tion that En{0) = 1 for all n. The simulation time for 
each initial condition was t = 5 x 10^ (~ 1.1 x 10“^ x Tq 
with Tq = V2tt). The results will be discussed below. 


A. Extreme Events: q 

Looking first of all at the distribution for the maxi¬ 
mum amplitudes of each extreme event (Fig. 5a), a clear 
trend emerges. Each distribution in k appears to exhibit, 
at least to a very good approximation, an exponential 
decay. That is, the probability of an extreme event at¬ 
taining ever higher amplitudes decreases exponentially, 
independent of the value of the coupling strength. In 
fact, after fitting these curves to determine the exact na¬ 
ture of the decay, it was established that the exponential 
function f{q) = aexp(—6g'^) is an appropriate fit, where 
the parameters a, 6, and c are constants. The results 
for this fit are contained in Fig. 6a. Evidently, the cou¬ 
pling strength n influences the value of the corresponding 
exponents: an increase in k sees also an increase in the 
exponent b. This manifests itself in the narrowing of 
the distributions seen in Eig. 5a. A simple explanation 
for this increase in b is that stronger couplings inhibit a 
unit’s ability to attain large amplitudes. In each case the 
parameter c remains close to c = 1, and even shows some 
convergence to this value. However, fitting the curves 
with c = I (not shown) produced significant errors, com¬ 
pared to the the case c ^ f, especially in the tails of the 
corresponding distributions. 

This fitting now enables some predictive power. It is 
possible now to predict the probability of an extreme 
event, with a certain amplitude (or greater than a given 
amplitude), taking place. In addition, one can compute 
the expected value of an extreme event E(qythres) = 
dq[q X aexp(—5(7'^)]. For k = 0.2,0.4, f.O one obtains 
E{qy,thres) = 5.06,4.69,4.33 respectively. 

The distributions for the duration of the extreme 
events (Fig. 5b) shows that events typically last no longer 
than t = 5 times units. The shape of this distribu¬ 
tion is maintained, regardless of the considered coupling 
strengths. However, the peak tends towards zero with in¬ 
creasing coupling strength. Moreover, the height of the 
peak increases, meaning that short duration events be¬ 
come more likely than the long duration events. From 
this result, it is not apparent that some of the localised 
structures have oscillatory behaviour, which can result in 
multiple extreme events in quick succession. One could 
argue that this succession of events instead accounts for 
a single extreme event, thus increasing the overall dura¬ 
tion of an event. However, this is extremely difficult to 
resolve in such a statistical analysis. 

Since there is no unique way of defining the wait¬ 
ing time between extreme events we consider two dif¬ 
ferent definitions which, interestingly, produce qualita¬ 
tively different results. The first considers extreme events 
when the consecutive events can happen anywhere on the 
chain. More precisely, if qn(ti) > qthres and qm(t2) > 
qthres where ti < 0 denote the start times of consecu¬ 
tive events, at chain locations n,m respectively, then the 
waiting time between events is simply truait = ^2 — ti- 

Secondly, we look at the waiting times between the 
end of one extreme event and the beginning of the next 
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extreme event that occur at the same chain location. To 
be more precise, if g„(ti) < qthres and g„(t 2 ) > qthres, 
where n denotes the chain location, ti marks the end of 
one event, and ^2 marks the beginning of the next, then 
the waiting time between these events is defined as above. 
Analogous definitions in terms of the energetic threshold 
can also be made. 

The distributions obtained when considering the first 
definition appear in Fig. 5c. They show a tendency for 
the waiting times to increase with increasing k. This 
is represented in the results by a lowering of the large 
peak and a broadening of the distribution. Interest¬ 
ingly though, additional peaks, with amplitudes increas¬ 
ing with K, emerge. To explain this, we should again refer 
to the heat maps contained in Fig. 1. For low k values 
there can simultaneously coexist multiple high amplitude 
structures, which are all likely candidates when search¬ 
ing for units that exhibit extreme events. In addition, 
the island structures consist of multiple units, and a unit 
crossing the threshold often pulls one or both of its neigh¬ 
bours beyond the threshold. Combined, these two factors 
result in the extreme events happening in quick succes¬ 
sion. Further, at low coupling strengths, it would seem 
that the oscillatory behaviour of the islands cannot ap¬ 
pear in the statistics because the frequency of oscillation 
is much slower than the rate at which new events occur. 
However, moving into the higher coupling regime, one 
sees that there are fewer island structures, and so less 
candidates for extreme events. This allows for the oscil¬ 
latory nature of the island structures to enter into the 
distributions. Of course, we still have the main peak due 
to the pulling effect mentioned above, but its amplitude 
is reduced because there are fewer island structures con¬ 
tributing to the short waiting times between consecutive 
events. 

When considering the second definition, the obtained 
distributions are even more complicated (see Fig. 5d). 
They exhibit a complex peak structure. Moreover, they 
seem to have heavier tails than the previously consid¬ 
ered distributions. In fact, this waiting time distribu¬ 
tion can be said to consist of two distinct regimes. The 
first regime is a correlated regime - the peak structure 
in Fig. 5d, which takes place on short time scales, and 
is almost certainly related to the ability of chain units 
to temporarily retain energy allowing for a successive se¬ 
quence of events. The second regime is an uncorrelated 
regime, taking place on long time scales (see Fig. 6-3), 
where new events happen at seemingly random intervals. 
This part of the distribution has a decaying exponential 
dependence. The emergence of the peak structure is dis¬ 
cussed in more detail in Sec. VI. 


B. Extreme Events: E 

The distributions of maxima in energy E for extreme 
events follow a similar behaviour to that of the ampli¬ 
tudes q in that they exhibit an exponential dependence 






Figure 5. Threshold crossing statistical distributions for q 
(left) and En (right). Row 1) maximum amplitudes for each 
event considered extreme; 2) event duration times for events 
considered extreme; 3) waiting times between events happen¬ 
ing anywhere on the chain; 4) waiting times between con¬ 
secutive events happening at the same chain location. More 
details in text. 


(see Fig. 6-2). However, in contrast, there is a broaden¬ 
ing of the distributions and an increase in the peak height 
with increasing k. This suggests that the chain preserves 
(maybe even enhances) its ability to localise energy when 
the strength of the interactions increase. These changes 
are not as pronounced as in the case of the q distribu¬ 
tions. 

The distributions for event duration also follow similar 
behaviour to their q counterpart, as the likelihood of long 
duration extreme events decreases with increasing k. In 
contrast though, a peak structure emerges. This is re¬ 
lated to the previously mentioned observation where the 
oscillations of the energy content of a unit which, while 
present, are often relatively small compared to oscilla¬ 
tions in q. This effect enters the statistics because of- 
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Figure 6. Color online. Lin-log plots for some of the results 
contained in Fig. 5 - panels 1),2),3),4) correspond respectively 
to panels a),e),d),h) from Fig. 5 (empty symbols - further de¬ 
tails are contained in the key). Also shown are fitted curves 
obtained using the function f{q) = aexp{—bq‘^) (filled sym¬ 
bols). The corresponding values b and c are contained in the 
key (c = 1 in the cases where c is not explicitly given in 
the key), while a is a constant representing the bin height at 
q = 4:, E = 3. The coupling strength related to each symbol 
is contained in the key in the top right panel. 


ten the energetic oscillations happen above the threshold 
- rather than crossing the threshold to end an extreme 
event, before recrossing it again to start another - thus 
increasing the duration of events. 

Looking at the time between events when using the 
first definition from the previous section (time between 
consecutive events occurring at any chain location), we 
see that the time decreases, with long waiting times be¬ 
coming less likely for larger coupling strengths. This dis¬ 
tribution is more regular than in the q threshold case, as 
there appears to be just a single maximum peak located 
at zero. 

However, these distributions become irregular when 
we consider the second definition for the waiting times 
between events (time between consecutive events occur¬ 
ring at each chain location), with additional peaks form¬ 
ing, and the tails becoming heavier. These complicated 
structures suggest that there are multiple mechanisms 
involved in the creation of high amplitude/energy fluctu¬ 
ations. 


VI. PEAK STRUCTURE: DISCRETE 
BREATHERS? 

Examining a sequence of extreme events (in q) one sees 
that the localised structures oscillate in time. This oscil¬ 
latory behaviour means that a single localised structure 
may cross the threshold qthres = 4 multiple times (see 
Fig. 7 - left panel) resulting in multiple extreme events 
coming from a single localised structure. 


It seems reasonable that this oscillatory behaviour may 
partially explain the interesting collection of peaks seen 
in the statistics of the previous section. Of course, the 
oscillatory behaviour is quite irregular. However, perhaps 
it is possible to predict the frequency of oscillation for 
these localised structures by exciting exact breathers of 
similar amplitude. 

It is worth remembering that the local Morse poten¬ 
tial is a soft anharmonic potential, which means that the 
larger the amplitude of the oscillation, the slower the fre¬ 
quency of oscillation. Therefore, the larger amplitude 
localised structures (that can be expected for low values 
of k) should oscillate with a lower frequency than the 
localised structures with lower amplitude. The same is 
true for the breather solutions. To obtain such breathers 
numerically one is referred to [20, 21]. 

Extracting the power spectrum from the system, for 
various values of the coupling strength k, one sees that 
a wide range of frequencies are excited. These frequen¬ 
cies include the windows of frequencies that would corre¬ 
spond to discrete breathers. However, the corresponding 
powers do not suggest that these windows are in anyway 
more significant than any other part of the spectrum. 
Therefore, one must conclude that discrete breathers are 
not playing a key role in the system dynamics. A better 
way to characterise the evolution of these high amplitude 
structures would be to characterise them as multi-site 
chaotic (both spatially and temporally) breathers. 

So if discrete breathers do not provide a solution to 
the emergence of the peaks that are present in the sta¬ 
tistical distributions of Sec. V, then where do they come 
from? To find an answer to this question, it is worth 
considering the evolution of a single initial condition at a 
time when one of the units undergoes multiple successive 
extreme events (see Fig. 7). One thing that is apparent 
in this figure is that prior to each extreme event there is 
a number k of subthreshold oscillations, i.e. oscillations 
whose amplitudes are less than qthres, with fe > 0. 

Taking a single initial condition and decomposing the 
waiting time statistics by the number of sub-threshold os¬ 
cillations between consecutive extreme events one finds 
that the obtained distributions indeed exhibit a peak 
structure: Fig. 7 - right panel. Moreover, this peak struc¬ 
ture (from a single initial condition) coincides extremely 
well with the peak structure in the distributions (from an 
ensemble of initial conditions) that are obtained when no 
such decomposition is carried out. 

Before summarising our findings it is worth mention¬ 
ing an aspect of our study that has gone largely unmen¬ 
tioned: the effects of chain length, relating to finite size 
effects, have not been considered, i.e how would chain 
length influence the statistical distributions? To answer 
this question, we computed the spatial autocorrelation 
function (see Fig. 8). The sharp decay in correlations 
for all chain lengths lead us to believe that finite size ef¬ 
fects are not an issue, and there would be no qualitative 
changes to any of the statistical distributions. 







Figure 7. Color online. Left: Time series for a single chain 
unit, which crosses the threshold gn = 4 multiple times. 
Right: The stacked histogram shows the distribution of wait¬ 
ing times between consecutive extreme events when decom¬ 
posed into the number of subthreshold oscillations between 
events. The red curve (solid line) shows the full distribu¬ 
tion (from previous section) when no such decomposition has 
taken place. This curve has been normalised such that its 
largest peak coincides with that of the stacked histogram. In 
both panels k = 0.2. 



Figure 8. Color online. Spatial autocorrelation function for 
amplitudes q (left) and energy E (right), with different chain 
lengths, and with coupling strength re = 0.1. 


high amplitude and the high energy fluctuations that can 
result because of the coupled dynamics. Some of these 
fluctuations cross a threshold thus instigating what are 
termed extreme events. With respect to the dynamics, 
we observed that energy could move along the chain and 
become localised allowing for the formation of high en¬ 
ergy island-like structures. Importantly, we have demon¬ 
strated that energy localisation should be considered in 
terms of a high energy background. These high energy 
island-like structures, possibly containing multiple units, 
can move coherently along the chain and influence the 
statistical properties of extreme events. 

A statistical analysis of the extreme events has unveiled 
some interesting properties of such events. For example, 
the distributions for amplitudes of extreme events ex¬ 
hibits a clear exponential dependence. In addition, the 
waiting times between consecutive events happening at 
the same chain location can be broken into regimes: the 
first a correlated regime where multiple events happen in 
quick succession due to a temporary retention of energy 
at some chain location, the second an uncorrelated regime 
where events happen seemingly at random. We have also 
demonstrated that the peaks in the distributions related 
to this first regime result because often there are a num¬ 
ber of subthreshold oscillations separating two extreme 
events. 

Certainly it is of interest to extend the present study to 
one where different underlying potentials are considered. 
For example, using a hard or harmonic potential instead 
of the soft potential used here will result in dynamics that 
are both quantitatively and qualitatively different [16, 
17]. Thus, we can expect that the statistics of extreme 
events will also be different. 


VII. CONCLUSIONS & DISCUSSION 
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